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SECTION  1 


INTRODUCTION 

In  the  investigation  of  crack  growth  behavior  using  an  interferometric 
displacement  gage  (IDG)  [1]  +  ,  an  empirical  relationship  between  crack  length 
and  measured  displacement  under  a  load  provides  the  basis  for  a  nonvisual 
method  of  determining  crack  length.  In  theory,  the  displacement  is  assumed  to 
be  a  homogeneous  linear  function  of  load.  Hence,  a  single  displacement 
measurement  would  be  sufficient  to  define  a  load-displacement  (P-6)  curve. 
However,  in  the  laboratory,  only  a  portion  of  the  load-displacement  curve  is 
linear. 

The  linear  portion  is  characterized  by  the  compliance,  C,  which  is 
the  inverse  of  the  slope,  shown  in  Figure  1.  This  figure  shows  a  typical 
load-displacement  curve  and  a  graphical  definition  of  the  compliance.  Instead 
of  using  compliance  which  has  units  of  length  per  force,  a  nondimens ional 
quantity,  EBC,  will  be  used  where  E  is  an  elastic  modulus  and  B  is  the 
thickness  of  the  specimen.  The  relationship  between  EBC  and  crack  length 
developed  in  this  report  can  be  applied  to  specimens  having  a  compact  type 
(CT)  configuration  and  to  materials  having  an  elastic  response.  In  the 
following  discussion,  expressions  are  developed  for  the  nondimensional 
compliance  or,  equivalently,  the  displacement  between  the  crack  surfaces  as 
a  function  of  distance  from  the  crack  tip  and  crack  length. 


Numbers  in  Brackets  Refer  to  References 


SECTION  2 


GENERAL  EXPRESSIONS  FOR 
REPRESENTING  DISPLACEMENTS 

Using  a  general  asymptotic  form  for  crack  surface  displacement  in  [2], 
the  nondimensional  compliance  can  be  written  as: 

EBC  =  A Vr  +  A  r  +  A  r3^2  +  ...  (1) 

U  x. 

where  r  is  the  distance  from  the  crack  tip  measures  along  the  crack  surface 
and  Aq,  A^,  ...  are  unknown  coef ficients.  The  number  of  coefficients  could 
be  terminated  arbitrarily  to  produce  any  suitable  degree  of  fit.  Also,  certain 
coefficients  could  be  chosen  to  have  specific  values.  From  theoretical 
considerations,  as  r  -*■  0,  the  first  coefficient  can  be  related  to  the  stress 
intensity  factor,  K, 


A 


0 


8B 

P/fn 


K, 


(2) 


where  B  and  P  are  the  thickness  and  applied  load,  respectively. 


K  =  /2ttx  a  (x,0), 
x-Kj  yy 


and  x  and  c  are  the  distance  and  normal  stress,  respectively,  ahead  of  the 

yy 


crack  tip. 
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Squaring  (1)  yields  another  form  of  a  displacement  expression 


(EBC) 2  =  AqT  +  A|r3/2  +  A^r2  +  ...  (3) 

where  A^,  as  in  (1),  and  A^,  A^»  ...  are  unknown  coeff icients.  Again,  the 
expression  can  be  terminated  at  any  degree  or  some  coefficients  can  be  defined 
a  priori. 

Two  representations  will  be  developed.  One  is  based  on  a  truncation 
of  the  series  expression.  The  other  is  based  on  a  modification  of  a  far  field 
representation  of  the  displacements.  The  truncated  series  representation  is 
explored  in  the  next  section. 


SECTION  3 


EVALUATION  OF  COEFFICIENTS 
FOR  TRUNCATED  SERIES 


In  the  development  of  an  equation  for  the  displacement,  now  expressed 
as  EBC,  for  a  wide  range  of  values  of  crack  length  and  distance  from  the  tip, 
Orange's  equation  [3]  which  was  developed  for  a  bend  type  specimen  was  examined. 
Orange  developed  an  equation  by  approximating  the  shape  of  the  crack  surface 
with  a  conic  section.  The  equation  has  the  form 


2 


=  _A_  (£)  +  _S!_ 

2+m  a'  2+m 


(4) 


where  a  is  the  crack  length,  6^  is  the  displacement  of  the  crack  surface  at 
r  =  a  under  a  given  load,  and  m  is  the  conic  section  coefficient.  Note  that 
only  two  nonzero  coefficients  in  (3)  are  used  to  model  (4).  From  analytical 
results  for  bend  and  single  edge  notch  (SEN)  specimens,  the  conic  section 
coefficient  can  be  expressed  as  a  function  of  crack  length  ratio 


m  =  -0.3  +  15  (a/W)n,  (5) 

where  n  is  2.3  for  bending  or  3.3  for  SEN,  and  W  is  the  width  of  the  specimen. 
However,  extending  an  expression  similar  to  (5)  for  CT  specimens  did  not 
provide  reasonable  results  for  the  displacement  or,  equivalently,  the 
compliance  lor  a  wide  range  of  values  of  a  and  r. 
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In  the  following  paragraphs,  the  coefficients  for  various  expressions 
of  compliance  will  be  evaluated  from  a  least  square  error  fit  to  analytical 
data  on  crack  opening  displacements.  A  finite  element  program  described  in 
[A]  was  used  to  generate  the  analytical  values  of  EBC.  The  results  are 
tabulated  in  Table  1  for  various  values  of  crack  lengths,  a/W,  and  r. 
Coordinates  and  notation  for  a  CT  specimen  are  shown  in  Figure  2. 

3.1  TWO  COEFFICIENT  FIT  FOR  (EBC)2 

Equation  (4)  is  a  specific  expression  which  incorporates  only  two 
unknowns,  6^  and  m,  of  the  more  general  expression  in  (3).  Even  though  m 
was  the  only  unknown  determined  by  an  empirical  fit,  both  5^  and  m  could  be 
determined  by  this  means.  Using  only  the  first  two  coefficients  of  integer 
powers  in  r  from  (3)  as  unknowns,  compliance  could  be  written  as 

(EBC)2  =  ar  +  gr2,  (6) 

where  a  and  3  are  unknown  coefficients.  Values  of  the  coefficients,  a  and  3, 
determined  from  the  fit  of  the  analytical  data  for  various  ratios  of  a/W  are 
given  in  Table  2. 

To  examine  the  validity  of  the  values  of  the  coefficients,  a  comparison 
between  (6)  and  (3)  was  made  for  r  ■*  0.  Using  (2),  one  obtains 
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TABLE  2 

VALUES  OF  a,  3,  and  (K’)2  VERSUS  a/W 


a/W 

0. 2500E+00 
0. 3000E+00 
0. 3500E+00 
0. 4000E+00 
0. 5000E+00 
0.6000E+00 
0. 7000E+00 
0. 8000E+00 


a 

0. 1525E+03 
0. 2017E+03 
0. 2598E+03 
0. 3343E+03 
0. 5724E+03 
0. 1102E+04 
0. 2678E+04 
0.9615E+04 


3 

0. 1437E+03 
0.2435E+03 
0. 3744E+03 
0. 5582E+03 
0. 1230E+04 
0. 2990E+04 
0. 9158E+04 
0 . 4954E+05 


(K 1 )  2 

0.1568E+03 
0. 2043E+03 
0. 2642E+03 
0. 3426E+03 
0. 6034E+03 
0.1206E+04 
0. 3004E+04 
0. 1098E+05 


For  a  CT  specimen,  the  stress  intensity  factor  is  given  by 


Bv'ft  W 


where  from  [5] 


f(S}  =  ~2J/3/2(-886+4-64($)~13‘32(S)2+14-72(f)3-5-6(S)4)' 

(1  W^ 


Defining  a  new  quantity. 


£<S>- 


and  using  (8) ,  (7)  can  be  restated  as 


a  =  A2  =  (K')2,  r  ■>  0. 


Since  a  was  determined  from  a  least  square  fit  to  numerical  results,  (11) 

may  not  be  satisfied  exactly.  Values  of  a  are  compared  with  values  of  (K') 

2 

in  Table  2.  When  a/W  is  small,  a  is  in  fairly  good  agreement  with  (K')  ,  b 

2 

as  a/W  increases,  values  of  a  deviate  more  from  (K.')  . 


Since  the  least  square  error  fit  using  (6)  follows  the  general  trend 

of  EBC ,  the  compliance  as  r  +  0  may  not  be  accurately  described.  The 

2 

deviation  in  a  from  (K')  occurs  because  (p)  is  forced  to  fit  larger  values 


of  EBC  away  from  the  crack  tip  as  a/W  increases.  Thus,  the  compliance 
behavior  further  away  from  the  crack  tip  becomes  more  dominant  in  the  fit 
and  the  resulting  equation  does  not  characterize  the  near  crack  tip  region. 


Additional  expressions  for  crack  surface  compliance  derived  from 
Equation  (1)  will  be  investigated.  In  the  following  equations,  the  near  tip 
behavior  is  enforced  analytically. 


3.2  TWO  COEFFICIENT  FIT  FOR  (EBC)2  BEYOND  CRACK  TIP  VICINITY 


From  the  theoretical  considerations  discussed  earlier,  the  first 
coefficient  in  (3)  is  assumed  to  be  given  by  (11).  Also,  it  was  decided 
that  'half-integer'  powers  should  be  considered.  To  maintain  simplicity, 
the  number  of  unknown  coefficients  should  be  kept  small.  To  generalize  (3), 
r  was  replaced  with  a  nondimensional  quantity,  p  =  r/W,  yielding: 


(EBC)2  =  (K')2  Wp  +  cqP3/2  +  c 1P2 


(12 


and 


(EBC)2  =  (K* ) 2  W„>  +  dQp2  +  dlP3 


(13 


where  c^,  c^,  d^,  and  d^  are  the  unknown  coefficients.  For  a  preliminary 
evaluation  of  the  coefficients,  the  data  for  a/W  =0.8  was  chosen  because  it 
appeared  to  be  the  most  difficult  to  fit.  The  results  of  fitting  (12)  and 


(13)  are  shown  in  Tables  3  and  4,  respectively.  The  comparison  of  the  root 
mean  square  deviation  for  the  two  fits  indicates  that  the  half-integer 
expression  provides  a  slightly  smaller  error.  Thus,  Equation  (12)  will  be 
investigated  further. 

3.3  TWO  COEFFICIENT  FIT  FOR  EBC  BEYOND  CRACK  TIP  VICINITY 

Since  a  direct  calculation  of  displacement  would  be  more  desirable 
than  a  squared  expression,  a  truncated  series  from  (1)  using  the  result  in 
(11)  and  expressed  in  terms  of  p  was  investigated, 

EBC  =  Ai  K'/p  +  bQp  +  bp'  (14 

where  b^  and  are  unknown  coefficients.  The  results  from  the  fit  using 

(14)  are  compared  to  the  results  using  (12)  in  Table  5.  Since  the  error 
using  (14)  was  generally  larger  than  the  error  using  (12),  the  coefficients 
in  (Id)  are  determined  for  other  values  of  a/W. 

3.4  COEFFICIENTS  AS  FUNCTIONS  OF  a/W 


Iron  (12) ,  values  of  c^  and  c^  were  determined  for  various  values  of 


a/W  and  art  tabulated  in  Table  6.  In  an  attempt  to  determine  convenient 


expressions  of  and  c^  as  a  function  of  a/W,  -Cq  and  c^  were  first  plotted 


a. ainsL  a/W  on  a  log-linear  scale.  It  appeared  that  -c^  and  could  vary 


with  respect  to  a/W  in  an  exponential  manner.  Polynomial  fits  of  the 


natural  log  of  -c^  and  e  as  a  function  of  a/W  were  developed.  Plots  of  the 


TABLE  3 

RESULTS  FROM  A  LEAST  SQUARE  ERROR  FIT 
FOR  a/W  =  0.8  USING  EQ.  (12), 

(EBC) 2=(K' ) 2Wp4C0p3/2+ClP2 


EBC 

EBC 

Percent 

r/W 

Anal . 

Calc . 

Error 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

0. 2502E-02 

0.6582E+01 

0 . 6566E+01 

-0. 2533E+00 

0.1000E-01 

0. 1341E+02 

0. 1335E+02 

-0. 4321E+00 

0. 1300E-01 

0. 1541E+02 

0 . 1534E+02 

-0.4198E+00 

0 . 1600E-01 

0. 1725E+02 

0. 1716E+02 

-0. 5237E+00 

0.2800E-01 

0.  2367E+02 

0 . 2345E+02 

-0. 9505E+00 

0. 4000E-01 

0. 2925E+02 

0. 2894E+02 

-0. 1055E+01 

0.7655E-01 

0. 4428E+02 

0.4381E+02 

-0. 1061E+01 

0.113IE+00 

0. 5798E+02 

0. 5761E+02 

-0. 6456E+00 

0. 1984E+00 

0. 8872E+02 

0. 8855E+02 

-0.1879E+00 

0. 2837E+00 

0. 1189E+03 

0. 1189E+03 

0. 6660E-01 

0. 4542E+00 

0. 1793E+03 

0.1794E+03 

0. 6460E-01 

0. 6248E+00 

0.  2398E+03 

0. 2397E+03 

-0. 1790E-01 

TABLE  4 


RESULTS  FROM  A  LEAST  SQUARE  ERROR  FIT 
FOR  a/W  =  0.8  USING  EQ.  (13), 

(EBC)2  =  (K,)2WfH-d0p2+d1p3 


r/W 

EBC 

Calc. -Anal . 

Anal . 

Calc . 

Anal 

(-) 

(-) 

(-) 

(%) 

0 . 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0. 2502E-02 

0. 6582E+01 

0 . 6629E+01 

0. 7125E+00 

0 . 1000E-01 

0 . 1341E+02 

0.1357E+02 

0 . 1182E+01 

0. 1300E-0I 

0. 1540E+02 

0.1561E+02 

0. 1320E+01 

0. 1600E-01 

0. 1725E+02 

0. 1747E+02 

0. 1299E+01 

0 . 2800E-01 

0. 2367E+02 

0.2391E+02 

0. 1009E+01 

0 . 4000E-01 

0 . 2925E+02 

0. 2951E+02 

0. 8845E+00 

0.7655E-01 

0 . 4428E+02 

0.4452E+02 

0 . 5402E+00 

0. 1131E+00 

0. 5798E+02 

0. 5830E+02 

0. 5553E+00 

0. 1984E+00 

0 . 8872E+02 

0.8897E+02 

0 . 2865E+00 

0. 2837E+00 

0. 1189E+03 

0. 1190E+03 

0 . 3970E+00 

0.4542E+00 

0. 1793E+03 

0. 1791E+03 

- . 9090E-01 

0 . 6248E+00 

0. 2398E+03 

0 . 2398E+03 

0 . 1670E-01 

COMPARISON  OK  RKSl'I.TS  FROM  LEAST  SQUARE  ERROR  FITS 
FOR  a / \<  =  0.8  USING  EQ .  (14), 
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TABLE  6 


VALUES  FOR  THE  COEFFICIENTS  RESULTING  FROM  EQ.  (12), 

2  2  3/2  2 

(EBC)  =(K')  Wp+cQ0  '  +c1P 


a/W 

fo 

fl 

0.30 

-24. 16573 

638.57910 

0.  35 

-38.79956 

981.62280 

0.40 

-69.99175 

1474.47400 

0.  50 

-231.55556 

3311.06909 

0.60 

-697.99036 

8111.06787 

0.70 

-1989.86206 

25256.97852 

0.80 

-7595.10645 

129157.88281 

I 


) 


data  and  results  of  several  polynomial  fits  to  the  data  are  shown  in 


Figures  3  and  4.  In  Figure  3,  the  quadratic  function  for  InC-c^)  vs  a/W 
appears  to  fit  the  data  much  better  than  the  linear  function.  In  Figure  4, 
the  quadratic  fit  for  In(c^)  vs  a/W  at  first  seemed  adequate,  but  after 
values  of  EBC  were  calculated  for  both  the  quadratic  and  cubic  functions,  the 
cubic  expression  was  chosen  due  to  the  significantly  improved  error. 

The  resulting  "best  fit"  expressions  for  c^  and  c^  vs  a/W  are: 

c0  -  exp[v0  +  vpf)  +  y,g>2]  (15) 


and 


c1  =  expl 


.■  (-) 

i  V 


+  vf)2  + 


'  3VW 


(§)3] 


(16) 


where  .  .  and  '  ,  are  constants. 

Oil  U  1  2  3 

3.6  EMPIRICAL  EQUATION  USING  TRUNCATED  SERIES 

In  summarizing  the  previous  results,  the  following  empirical  equation 

1/2 

wi.ic:.  was  based  on  a  truncated  series  in  .  was  found  to  best  represent  the 
displacement  along  th.e  crack  surface  as  a  function  of  crack  length  for  a  CT 
spoc imen : 


t 


Compa r i son  of  Linear  and  Quadratic  Fits  < 


(17 


EEC  =  [(^-f(J))Vg(f)p3/2+h(§)p2]1/2 

/2tt 

where  p  =  r/W,  nondimensional  distance  behind  the  crack  tip, 

f(§)  =  - ^3/2(.886+4.64(J)-13.32(J)2+14.72(f)3-5.6(g)4),  (9) 

(i^p 

g(f)  =  -exp[v0  +  Vi(f)  +v2(|)2],  (18 

h(f)  =  exP[60  +  63(g)  +  62(|)  +  63(g)3],  (19 


vQ  =  0.0285773,  6Q  =  1.79452, 

V1  =  10.0701,  63  =  24.1346, 

V2  =  1.25919,  &2  =  -37.1842, 

63  =  28.2394. 

A  comparison  of  the  desired  values  of  EBC  and  the  calculated  results  using 
Equation  (17)  is  shown  in  Table  7.  From  the  results  in  Table  7,  the  error  is 
no  longer  greater  than  ±  3%  over  the  ranges  0  _<  p  <  a/W  -  0.175  and 
0.25  <  a/W  <  0.80. 


TABLE  7 

RESULTS  FROM  THE  TRUNCATED  SERIES  EXPRESSION, 
EQ.  (17),  TO  CALCULATE  EBC 


/W  =  0.25 
a=r/W 


Anal. 

(-) 


Calc . 
(-) 


Calc . -Anal 
Anal . 
(%) 


.000000 

0.000000 

0.000000 

- 

.000781 

0.429670 

0.432968 

0.7676 

.003130 

0.858010 

0.867738 

1 .1338 

.004067 

0 . 982720 

0.989578 

0.6979 

.  005003 

1 .099800 

1.098220 

-0.1437 

.008756 

1 . 435200 

1 .456174 

1.4614 

.012508 

1.735400 

1 .744700 

0.5359 

.016413 

1.997800 

2.003759 

0.2983 

.020311 

2.229500 

2.234934 

0.2437 

.029416 

2.720000 

2.706385 

-0.5005 

.038514 

3. 129800 

3.116203 

-0.4344 

.056717 

3.859000 

3.829089 

-0.7751 

.074921 

4.463100 

4.455533 

-0.1695 

/W  =  0.30 


P=r/W 


Anal . 
(-) 


Calc. 

(-) 


Calc . -Anal 
Anal . 
(%) 


.000000 

0.000000 

0.000000 

- 

.  000940 

0.541600 

0.542119 

0.0959 

.003752 

1.082200 

1 ,085152 

0.2728 

.  004879 

1 . 238800 

1.238431 

-0.0298 

.  006006 

1 . 384600 

1 .375191 

-0.6795 

.010514 

1 .815300 

1 . 826007 

0.5898 

.015016 

2 . 196400 

2 .190332 

-0.2763 

.  021898 

2.680700 

2.660598 

-0.7499 

.02877f, 

3 . 089600 

3 . 067915 

-0,7019 

.044825 

3.933100 

3.882429 

-1 .2883 

.060876 

4.643900 

4 ,586820 

-1 .2291 

.092952 

5.915300 

5 . 821053 

-1 .5933 

.  125079 

7,004700 

6.927941 

-1 .0958 

TABLE  7  (Continued) 

RESULTS  FROM  THE  TRUNCATED  SERIES  EXPRESSION, 
EQ.  (17),  TO  CALCULATE  EBC 


a/W  =  0.35 


p=r/W 

EBC 

Calc. -Anal . 

Anal . 

Calc. 

Anal . 

(-) 

(-) 

(-) 

(%) 

.000000 

0.000000 

0.000000 

_ 

.001092 

0.664710 

0 . 664469 

-0.0363 

.004375 

1 . 332700 

1  .  333024 

0.0243 

.005686 

1 . 524200 

1 .521432 

-0 .1816 

.006997 

1 .701400 

1  .  689752 

-0.6846 

.012248 

2.238700 

2.246811 

0.3623 

.017492 

2.708900 

2 .699108 

-0.3615 

.027314 

3.439700 

3 . 406576 

-0.9630 

.037130 

4.047100 

4.011774 

-0.8729 

.  060044 

5.283500 

5.221133 

-1 .1804 

.  082984 

6.342500 

6 . 277997 

-1.0170 

.  128762 

8.273600 

8 .  161949 

-1 .3495 

.174603 

9 . 979900 

9.891178 

-0.8890 

a/W  =  0.40 


.  =r/W 

EBC 

Calc. -Anal. 

Anal  . 

Calc . 

Anal . 

(-) 

(-) 

(-) 

(%) 

. 000000 

0 . 000000 

0 ,000000 

.  o'- 

. 001251 

0 ,812070 

0.809429 

-0.3253 

. 005003 

1 . 626000 

1 .623625 

-0,1461 

. 006502 

1 . 859800 

1 . 853545 

-0.3363 

. 008000 

2 . 074700 

2.059201 

-0 ,7471 

.014000 

2 . 739200 

2.741730 

0.0924 

.  020000 

3.314400 

3 . 299209 

-0 .4583 

.032800 

4 . 333900 

4 . 287704 

-1 .0659 

.045594 

5 . 179900 

5 . 130266 

-0.9582 

.075429 

6 . 899500 

6.824113 

-1 .0926 

.  105333 

8 .398800 

8,327766 

-0 . 8458 

.  165016 

1 1 .184000 

11 .061826 

-1 .0924 

.224762 

13.715000 

1 3 . 6231 98 

-0.6694 
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TABLE  7  (Continued) 

RESULTS  FROM  THE  TRUNCATED  SERIES  EXPRESSION, 
EQ.  (17),  TO  CALCULATE  EBC 


a/W  =  0.50 


p= r/W 

EBC 

Anal. 

Calc. 

Calc. -Anal. 
Anal . 

(-) 

(-) 

(-) 

(%) 

.  000000 

0.000000 

0.000000 

- 

.001562 

1 .206500 

1 .198247 

-0.6840 

. 006254 

2.425100 

2.406589 

-0.7633 

.008133 

2.772800 

2.750034 

-0.8210 

.010006 

3.091400 

3.056865 

-1 . 1171 

.017511 

4.108900 

4.081822 

-0.6590 

.025016 

4.975400 

4.927234 

-0.9681 

.043771 

6 . 780500 

6 . 684945 

-1 .4093 

. 062521 

8.293300 

8.192843 

-1.2113 

.106286 

11 .419000 

1 1 .296827 

-1 .0699 

. 150032 

14.245000 

14.130841 

-0.8014 

. 237587 

19.663 000 

19.475500 

-0.9536 

.325079 

24.812000 

24.626057 

-0.7494 

a/W  =  0.60 

p=r/W 

EBC 

Calc. -Anal. 

Anal . 

Calc. 

Anal . 

(-) 

(-) 

(-) 

(%) 

.000000 

0 . 000000 

0 . 000000 

- 

.001873 

1 .873700 

1.850269 

-1 .2505 

. 007498 

3.776300 

3.717831 

-1 .5483 

.009752 

4 .318900 

4.251172 

-1 .5682 

.012000 

4.815500 

4.729215 

-1 .7918 

. 021003 

6 . 446100 

6.336740 

- 1 .6965 

. 030000 

7.825000 

7.677175 

-1 .8891 

,054673 

1 1 . 000000 

10,768578 

-2 . 1038 

. 079365 

13,706000 

13.470174 

-1.7206 

. 136889 

19.431999 

19. 178482 

-1 .3046 

.194476 

24 . 784000 

24 . 553379 

-0 , 9305 

. 309651 

35.285999 

34 .946835 

-0.9612 

, 424762 

45 .576000 

45 . 1 60328 

-0 .9120 

TABLE  7  (Concluded) 

RESULTS  FROM  THE  TRUNCATED  SERIES  EXPRESSION, 
EQ.  (17),  TO  CALCULATE  EBC 


/W  =  0.70 


p=r/W 

EBC 

Calc.  -An. 

Anal. 

Calc. 

Anal 

(-) 

(-) 

(-) 

(%) 

.  000000 

0.000000 

0 . 000000 

.  002190 

3.221000 

3 .154046 

-2.0787 

.008756 

6.507900 

6.351496 

-2.4033 

.011384 

7.451700 

7.272814 

-2.4006 

.014006 

8.316400 

8.103281 

-2.5626 

.024514 

11 .237000 

10.931138 

-2.7219 

.  035016 

13,717000 

13.333463 

-2.7961 

.065651 

19.899000 

19.366501 

-2 . 6760 

.  096254 

25.305000 

24.783155 

-2.0622 

.  167746 

37.063999 

36.634975 

-1  .  1575 

.239238 

48.354000 

48 . 083637 

-0.5591 

.  382159 

70 . 839996 

70 . 633438 

-0.2916 

.525079 

93 . 247002 

93,067078 

-0.1930 

/W  =0.80 


.  =r/W 

EBC 

Calc. -An 

Anal . 

Calc. 

Anal 

(-) 

(-) 

(-) 

(%) 

000000 

0.000000 

0 . 000000 

_ 

002502 

6 . 582500 

6 . 468560 

-1 .7309 

010000 

13 . 408000 

13. 162839 

-1 .8285 

013003 

15.405000 

15. 132161 

~1 .7711 

016000 

17 . 246000 

16.925253 

-1 .8598 

028000 

23.672001 

23 . 154379 

-2 . 1845 

040000 

29 .247999 

28.601858 

-2 . 2092 

076571 

44 . 282001 

43.392429 

-2 . 0089 

1  13079 

57,981998 

57 , 1  12049 

-1 ,5004 

198349 

88,718002 

87 . 953712 

-0.8615 

283683 

118.870003 

1 18.289337 

-0 .4885 

454222 

179 . 250000 

178.560394 

-0 . 3847 

624762 

239 . 759995 

238 . 760498 

-0.4169 

SECTION  4 


EVALUATION  OF  A  LINEAR  EXPRESSION 
WITH  CORRECTION  TERM 


An  expression  for  Che  crack  surface  displacement  as  a  function  of 
crack  length  for  a  CT  specimen  has  been  developed  by  Saxena  and  Hudak  [6]. 
The  form  of  the  expression  is  desirable  because  of  its  simplicity.  In  terms 
of  compliance,  the  expression  was 


EBC 


XR/W  -  X/W 
XR/W  +  0.25  EBC0 


(20) 


where  X  and  X^  are  the  distances  from  the  load  line  to  the  location  of 
displacement  on  the  crack  surface  and  to  the  axis  of  rotation*,  respectively, 
and  CQ  is  the  compliance  at  the  crack  mouth.  In  [6],  values  of  X^  are 
tabulated  for  various  values  of  a/W  and  C^  is  given  by  the  following  empirical 
f unction 


EBC  = ( 1+— ) (~ ^tS)2[1. 61369+12. 6778(a/W)-14 . 2311 (a/W) 2 
0  a/W  1-a/W 


(21) 


-16.6102(a/W)3+35.0499(a/W)4-14.4943(a/W)5] . 


For  a  given  a/W,  (20)  is  linear  in  X  or  in  r  since  X  =  a  -  r. 


The  axis  of  rotation  is  defined  in  [6]. 


Essentially,  (20)  is  a  good  representation  of  the  crack  surface 
displacements  in  the  vicinity  of  the  load  line.  A  plot  of  the  finite 


element  data  in  Figure  5  for  a  relatively  short  crack,  a/W  =  0.25,  in  a  CT 
specimen  shows  that  the  linear  behavior  is  a  good  approximation  to  the 
compliance  in  the  vicinity  of  the  load  line,  i.e.,  0.1  <  p  <  0.6.  However, 
note  that  the  linear  aporoximation  does  not  pass  through  the  crack  tip  where 
the  compliance  should  be  zero. 


4.1  CORRECTION  TERM 


The  addition  of  a  correction  term  to  the  linear  compliance  expression 
would  be  necessary  to  obtain  the  near  tip  behavior.  Thus,  a  representation 
for  a  general  compliance  expression,  cf.  Equation  (1),  could  be  written  as 


C/CQ  =  aQ  +  a^  p  -  (correction  term), 


where  a  and  a^  are  functions  of  a/W  and  the  correction  term  is  a  function  of 

a/W  and  .  . 


Comparing  the  linear  terms  of  (20)  when  X  =  a  -  r  and  (22),  the 
coefficients  in  (22)  could  be  approximated  by 


XR/W  -  a/W 


a0  X^/W  4-0.25  and  al  X^/W  +  0.25  ' 


Values  of  a^  and  which  are  calculated  from  (23)  using  the  values  of  X^ 


.’hich  are  given  in  [6]  are  listed  in  Table  8. 
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TABLE  8 


ESTIMATED  VALUES  OF  THE  COEFFICIENTS  FOR 
THE  LINEAR  TERMS  IN  EQ.  (22) 


a/W  aO  ai 

0. 2500E+00  0.1380E+00  0.1720E+01 
0. 3000E+00  0. 2040E+00  0.1590E+01 
0. 3500E+00  0.1070E+00  0.1490E+01 
0. 4000E+00  0.9190E-01  0.1400E+01 
0. 5000E+00  0. 6320E-01  0.1250E+01 
0. 6000E+00  0. 4330E-0I  0.1130E+01 
0. 7000E+00  0. 3710E-01  0.1014E+0I 
0. 8000E+00  0. 2960E-0I  0.9240E+00 


In  developing  the  correction  term,  the  theoretical  behavior  of  the 


compliance  as  p  -*•  0  which  is  described  in  (2)  and  (1)  will  be  imposed.  From 
the  result  in  (2)  and  using  (8)  and  (10) ,  the  correction  term  in  (22)  must 
have  the  following  behavior. 


(correction  term)  -  a  =  - 


K*  /W 
EBC„ 


0. 


(24) 


The  correction  term  must  also  approach  zero  far  from  the  crack  tip  so  that 
the  form  of  the  Saxena-Hudak  expression  is  recovered.  Various  selections  for 
the  functional  form  of  the  correction  term  are  discussed  later. 


4.2  SOLUTION  TECHNIQUE 

The  procedures  which  are  used  to  determine  the  unknown  quantities  in 
(22)  involve  a  fit  to  the  known  data  and  an  application  of  a  measure  of  the 
nearness  of  the  fit.  Since  (22)  is  expected  to  be  nonlinear  with  respect  to 
the  unknown  coefficients,  a  direct  solution  is  not  generally  possible.  Two 
types  of  iterative  techniques  are  considered  to  obtain  a  solution  for  the 
unknowns . 


One  technique  is  to  select  initial  values  for  those  unknown  quantities 
which  contribute  the  nonlinear  terms  in  (22)  and  solve  in  a  least  square  error 
sense  the  unknowns  in  the  linear  terms.  Then,  update  the  values  of  the 
coefficients  in  the  nonlinear  terms  with  the  most  recent  solution  and  repeat 
the  least  square  error  solution.  The  iteration  was  terminated  when  the 
solution  converges  to  within  a  small  fixed  percent  of  the  preceding  solution. 
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A  second  technique  which  is  basically  a  trial  and  error  method  uses 
an  initial  guess  for  the  solution  of  the  unknowns  and  then  the  deviations  and 
percent  errors  in  the  calculated  compliances  are  evaluated.  The  initial  guess 
is  incremented  and  a  new  set  of  deviations  are  calculated.  Successive  incre¬ 
ments  are  made  to  the  values  of  the  unknowns  until  the  deviations  are 
acceptable . 

Two  criteria  are  applied  to  the  deviations.  For  one  criterion,  the 
root  mean  square  deviation  (RMSD)  is  examined  for  a  minimum.  For  the  second 
criterion,  the  distribution  of  the  percent  errors  of  the  deviations  is 
examined  for  uniform  scatter  about  the  desired  values.  Generally,  these  two 
criteria  cannot  be  satisfied  simultaneously.  However,  a  solution  for  the 
unknowns  is  selected  when  the  RMSD  is  near  a  minimum  value  in  order  that  the 
maximum  error  in  the  scatter  is  reduced  to  a  small  value.  A  result  based  on 
this  subjective  technique  can  be  obtained  rather  quickly  if  the  initial  guess 
is  close  to  a  minimum  RMSD.  An  example  of  the  RMSD  using  this  procedure  is 
given  in  Section  4.4. 

A. 3  FINITE  ELEMENT  RESULTS  FOR  COMPLIANCE 

The  unknowns  in  (22)  are  determined  so  that  the  calculated  values  of 
compliance  from  (22)  agree  with  known  values  of  compliance.  Since  Newman's 
results  [7]  have  been  used  to  develop  various  expressions  for  compliance  or 
displacement  in  the  vicinity  of  the  load  line  for  a  CT  specimen,  those  results 
•  ire  ..•■•mb  ined  with  the  finite  element  results  determined  for  this  study  for  the 
intermediate  distances  between  the  load  line  and  the  crack  tip.  These 
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combined  values  are  shown  in  Table  9.  It  is  noted  that  the  finite  element 
results  from  the  analysis  used  for  this  report  did  not  significantly  differ 
from  Newman's  results  where  the  two  sets  of  results  overlapped.  The  data 
listed  in  Table  9  are  used  to  determine  the  coefficients  in  the  subsequent 
expressions  for  compliance. 


4.4  EXPRESSIONS  FOR  THE  CORRECTION  TERM 


Various  functional  expressions  for  the  correction  term  in  (22)  were 
evaluated.  For  all  the  expressions,  the  near  crack  tip  behavior  given  in  (24) 
is  satisfied  analytically.  One  of  the  expressions  is 


c  ao  ,  ao 

=  V  +  V  "  T  6 


where  \  =  K'v'W/(EBC^)  and  n  is  an  unknown  constant  which  was  selected  to 
have  integer  values.  The  first  solution  technique,  described  in  Section  4.2, 
was  used  to  determine  the  coefficients,  a^  and  a^,  with  their  initial  values 
given  in  Table  8.  Solutions  were  obtained  but  the  values  of  a^  and  a^ 
differed  significantly  from  the  initial  values  such  that  if  values  of  p  were 
larger  than  those  used  in  region  of  fit,  the  solution  was  not  realistic. 


Another  expression  for  the  correction  term  in  (22)  yields 


— -  =  a  +  a  o  -  f  (A)/g(>  p)  , 
L0  U 


where  f  and  g  are  polynomial  functions  in  /p.  Only  polynomial  equations 
involving  no  more  than  two  coefficients  in  each  function  were  investigated. 
Again,  the  first  solution  technique  in  Section  4.2  was  used  to  determine  the 
unknown  coefficients  in  (26).  Generally,  the  results  from  the  least  square 
error  solution  did  not  fit  the  intermediate  compliance  data  very  well.  Also, 
the  correction  term  did  not  approach  zero  rapidly  enough  as  0  increased. 

Since  the  exponential  behavior  appeared  to  produce  a  better  fit, 
various  other  expressions  involving  exponential  behavior  were  tried.  The 
following  expression  for  the  correction  term  produced  a  good  fit 


(correction  term) 


(l+a2)aQ 


1+a 


1+a2  exp[a„a 


2A/p  ] 


2  0 


(27) 


where  an  is  also  an  unknown  quantity.  For  any  value  of  a2>  the  near  crack 
tip  behavior  in  (24)  is  satisfied.  To  minimize  the  number  of  coefficients 
which  arc  functions  of  a/W  in  (27),  i.e.,  to  maintain  a  simple  expression  for 
the  compli.ince,  is  taken  to  be  a  constant. 

Using  the  data  for  a/W  =  0.5,  several  values  of  were  chosen  and 
trie  second  solution  technique  in  Section  4.2  was  used  to  determine  a^  and  a^. 
The  second  technique  was  applied  because  the  first  technique  did  not  yield 
acceptable  results  for  the  exponential  expression  in  (27).  The  value  of 

The  value,  a0  = 

i. 


which  produced  a  good  result  was  approximately  2. 


2 ,  was 


As  an  example  of  the  application  of  the  second  solution  technique 
when  (27)  is  used  in  (22) ,  increments  were  subtracted  and  added  to  the 
initial  values  of  a^  and  a^  which  were  chosen  from  Table  8  for  a  given  a/W 
value.  Compliance  values,  percent  errors  and  a  RMSD  were  calculated  for 
each  combination  of  and  a  .  An  array  of  RMSD  was  created  for  the 
values  of  and  the  values  of  a^  as  illustrated  in  Table  10.  The  percent 
errors  in  calculated  compliances  were  examined  for  the  values  of  a^  and  a^ 
in  the  vicinity  of  the  minimum  RMSD.  If  a  set  of  percent  errors  was  uniformly 
distributed  within  a  ±  3%  scatter  band  then  the  values  of  a^  and  for  that 
set  were  selected  as  a  solution. 

If  none  of  the  scatter  bands  were  within  ±  3%,  then  smaller  increments 
would  be  applied  to  the  values  of  a^  and  a^  at  the  smallest  RMSD  value  to 
define  more  accurately  the  values  of  a^  and  a^  for  a  minimum  RMSD.  Again, 
the  percent  errors  would  be  examined  for  a  uniform  distribution  of  errors 
within  a  ±  3%  scatter  band.  The  procedure  was  repeated  until  a  ±  3%  scatter 
band  was  achieved  or  the  smallest  increments  of  change  were  0.0005  for  a^  and 
0.005  for  a,.  Then,  values  of  a.  and  a,  were  selected.  The  values  of  a_  and 

i  0  1  0 

a^  which  were  chosen  as  solutions  for  each  a/W  are  listed  in  Table  11. 

Next,  Uq  and  a^  were  considered  as  functions  of  a/W.  A  least  square 
error  polynomial  fit  was  used.  Figures  6  and  7  show  a^  and  a^  plotted 
against  a/W  to  illustrate  the  behavior  of  the  functions.  A  second  order 
polynomial  proved  satisfactory  for  a^  and  a  third  order  polynomial  was  needed 
for  a^ .  The  order  of  the  polynomials  was  kept  small  as  possible  to  preserve 
the  simplicity  and  yet  maintain  the  accuracy  of  the  fit. 
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TABLE  10 

ARRAY  USED  TO  EVALUATE  a0  AND  aj  FOR  EQ .  (27) 
T1TUTEI)  INTO  (22)  WHEN  a/W  =  0.4 


TABLE  11 


FINAL  VALUES  OF  THE  COEFFICIENTS  FOR  THE 
LINEAR  TYPE  EQUATION  WITH  CORRECTION  TERM 


a/W 

0. 2500E+00 
0. 3000E+00 
0. 3500E+00 
0. 4000E+00 
0. 5000E+00 
0.6000E+00 
0. 7000E+00 
0. 8000E+00 


0.8800E-01 
0.9320E-01 
0.9450E-01 
0. 9250E-01 
0 . 8250E-01 
0.6900E-01 
0. 6100E-01 
0. 4230E-01 


0.1820E+01 
0. 1637E+01 
0. 1485E+01 
0. 1375E+01 
0. 1205E+01 
0. 1090E+01 
0. 9750E+00 
0. 9085E+00 


I 


u  I 


4.5  EMPIRICAL  EQUATION  USING  CORRECTION  TERM 


The  relationship  developed  in  this  section  has  the  final  form: 


C  _ 

c  an  ai  p  —  _ 

0  l+2*exp[1.5A/p/aQ] 


where 


Cg  is  the  compliance  at  the  crack  mouth  opening,  Eq.  (21), 
a  =  K’/W/(EBC0) , 

aQ  =  .024  +  . 457 (a/W)  -  .897(a/W)2  +  .445(a/W)3, 
a_L  =  2.696  -  4. 279(a/W)  +  2.582(a/W)2, 


K'  is  given  by  (10)  and  (9). 


A  comparison  of  the  desired  values  of  EBC  and  the  calculated  results 
using  (28)  is  shown  in  Table  12.  From  the  results  in  Table  12,  the  error  is 
no  greater  than  ±  4%  over  the  ranges  0  <_  p  <_  a/W  +  0.25  and  0.25  a/W  _<  0.80. 
Note  that  (28)  produces  a  slightly  larger  deviation  from  the  desired  values 
than  (17)  but  (28)  covers  a  wider  range  of  distances  from  the  crack  tip  than 
(17).  Also,  (28)  is  expected  to  be  applicable  for  values  of  p  beyond 
a/W  +  0.25. 


I  A 


i  i 


i-, 


Si  .  J 


I...J 


r 
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TABLE  12 


RESULTS  FROM  THE  LINEAR  EXPRESSION 
PLUS  CORRECTION  TERM,  EQ .  (28),  TO  CALCULATE  EBC 


a/W  =  0.25 


p=r/W 

EBC 

Calc . -Anal . 

Anal . 

Calc. 

Anal . 

(-) 

(-) 

(- 

) 

(%) 

.  .  .  \ 

IJ  «.  •  *.• 
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00  00  0 
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1  3  7  6  r- 
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0*i 

o 

.)  t 

8  7  3  Z  7 

J.  .  '  7  '  /  O  _? 

.  ,  (.  0  i  • 

.  7  t>  -  /  A- 

•:  , 

(;>  ?  ’  j  7  3 

1  ,2240 ... 

■  /  '>_■  j 

; ,  ?  ?  8  o 

;  ♦ 

i  0  2  1  5 

•  J  *  2  |  _.i  9  H 

■  «  mj  i  2  3 

"  4  1  -  • 

J.  < 

1 4  e c ' 

0 , 9 ",  - 1  • 

■  ,  0  3.  0  3 

.w  »  v.  ...'.  7  ...:  V 
^  v  /  j  r.  i’ )  f  . 

07  7  8  o 

.  1-!  2  v" 

■  .  2  0  0 

3  t  8  3  ?  0  0 

.  !  , 

20  418 

0  v  7  i  .2 

'!  ••  V  .;  v  0 

.  \*  - 

21  iOi 

-  1  .  -  C  7  2  f: 

RMSD 

ERROR  =  0.1066 

a/W  =  0. 30 

j=r/W 

(-) 

EBC 

Anal.  Calc. 

(-)  (-) 

Calc. -Anal. 
Anal . 

(%) 

-  r.  • ,  - 

■■  '(''CO  0,00000 

*".  A  :  .*■  •:>  0  ,  '.VO  7  77 

1  -  °  .*■  '  " 

».  .  ..  •  •  • 

■  -  74 J  : '  S;7 ' 

'\  :  * 

-  .  -  . :•  \  .  ’.  8  0  4  1 

•  0  •  3.  i  • 

_  f  V  ’ "  *  . 

< 

-  -  7  ’  f  -  1  , 

-  :i  8  1  ■ 

RMSD  ERROR  =  0.0995 
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TABLE  12  (Continued) 


RESULTS  FROM  THE  LINEAR  EXPRESSION 
PLUS  CORRECTION  TERM,  EQ.  (28),  TO  CALCULATE  EBC 
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EBC  Calc. -Anal 


Anal . 

Calc. 

Anal . 

(-) 

(-) 

(%) 

.  • 

:■ .  ooooo 

-  o  6  •'}  7  1 

V,c:.736t 

■  •  77s  15 
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0  :  4  0 
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RMSD  ERROR  =  0.1703 
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EBC 

Calc. -Anal 

Anal . 

Calc. 

Anal . 

(-) 

(-) 

(%) 

RMSD  ERROR  =  0. 2237 


TABLE  12  (Continued) 

RESULTS  FROM  THE  LINEAR  EXPRESSION 
PLUS  CORRECTION  TERM,  EQ,  (28),  TO  CALCULATE  EBC 

a/W  =  O.  50 

p=r/W  EBC  Calc. -Anal. 

Anal.  Calc.  Anal. 

(-)  (-) _  M _  (%) 

■ .  :•  -  ,■  o  ~  "  ,  ■ ;  v  o  o .  o  o  o  o  o 


RMSD  ERROR  =  0.3131 


a/W  =  0.60 


»-=r/W 
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Anal. 

Calc. 
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RMSD  ERROR  =  0.7538 


TABLE  12  (Concluded) 


RESULTS  FROM  THE  LINEAR  EXPRESSION 
PLUS  CORRECTION  TERM,  EQ.  (28),  TO  CALCULATE  EBC 


a/W  =  0.70 


<*1 


v  =r/W 


Anal. 

(-) 


Calc. 

(-) 


Calc . -Anal . 
Anal. 
(%) 


260  Of' 


'i  b  0  8  o 


1 . 2  i  3  . 


•J  ,  */  ...  r 


RMSD  ERROR  =  1.3005 


a/W  =  0.80 


EBC 

Anal.  Calc. 

(-)  (-) 


Calc. -Anal . 
Anal . 
(%) 


L.’J 


RMSD  ERROR  =  3.2479 


SECTION  5 


APPLICATION 


Two  equations  have  been  developed  for  compliance  as  a  function  of 
crack  length  for  a  CT  specimen.  For  indirect  crack  length  determination,  it 
is  desired  to  have  crack  length  expressed  in  terms  of  compliance  since 
compliance  is  measured  in  the  experiment.  A  direct  mathematical  inversion 
is  not  easily  obtained  because  both  expressions  are  complicated  functions  of 
i.  A  curve  fit  of  crack  length  as  a  function  of  both  compliance  and  measure¬ 
ment  location  could  be  tedious  and  difficult.  Thus,  a  numerical  method, 
namely  Newton's  Method,  is  applied  to  Equation  (28)  to  calculate  values  of 
crack  length  when  compliance  is  known. 

The  problem  is  to  find  a  when  a  compliance  value  C*  is  measured  at  a 
known  Location  d,  see  Figure  2.  From  the  form  of  (28),  one  can  write 

C*  =  C(a,  a  -  6), 

where  .  =  •..»  -  £  and  £  =  d/W,  or 


C*  =  C(a) , 


(29) 


where  C  represents  a  function  of  only  crack  length.  An  approximate  solution 
for  i  in  (29)  can  be  found  using  Newton's  Method.  An  iterative  solution  of 


(29)  becomes 


a. 

1 


c*-c(oti_l) 


i-1,  2,  ... 


(30) 


i-1 


d  /s 

— tc] 

da  a=a 


i-1 


Since  C(0  is  a  complicated  expression,  the  derivative  was  estimated  in  the 
numerical  calculations  as  follows 


d  r;-,  C(a  +  A/2)  -  C(a  -  A/2) 
dx[  J  %  A 


(31) 


where  *  0.05  rt. 

To  find  a  from  (30)  when  (31)  is  substituted  and  the  compliance  has 
been  determined,  an  initial  estimate,  a  =  is  needed.  Then,  the  iteration 

is  terminated  when  |a  -  a_.  is  sufficiently  small. 

in  the  beginning  of  a  crack  growth  test,  a  could  be  determined  from 
optical  measurements.  As  the  crack  grows,  a  value  of  the  initial  crack  length 
estimate  for  a  subsequent  compliance  value  can  be  the  previously  determined 
crack  length. 

The  results  of  the  iterative  procedure,  using  (30)  with  the  approxi¬ 
mation  for  the  derivative  given  by  (31),  are  shown  in  Table  13.  For  these 
examples,  values  of  compliance  were  chosen  from  the  analytical  values  for 
i  =0.8  used  to  obtain  the  compliance  express: jn.  In  all  examples,  the 
initial  value  for  ot^  was  0.7.  Note  that  the  number  of  iterations  for 
convergence  to  a  solution  lx*  is  relatively  small.  The  numerical  iteration 
procedure  was  terminated  when  ja_^  -  m  ^  |  <  0.0001. 
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SECTION  6 


SUMMARY 

Two  empirical  equations  were  developed  for  the  displacement  between 
the  crack  surfaces  of  a  CT  specimen.  The  equations  were  derived  from  curve 
fits  to  finite  element  results  which  were  determined  from  a  linear  elastic 
analysis.  Since  a  linear  analysis  was  used,  the  displacements  for  any  load 
were  presented  using  compliance,  the  ratio  of  the  change  in  displacement  to 
the  change  in  load.  To  avoid  extensive  curve  fitting  and  to  maintain  some 
simplicity,  the  number  of  unknown  coefficients  used  to  fit  the  data  for  a 
given  crack,  length  was  limited  to  two.  Both  equations  yielded  the  correct 
analytical  behavior  as  the  crack  tip  was  approached. 

One  equation  was  developed  from  a  truncated  series  expansion  for 
the  square  of  the  compliance  in  terms  of  integer  powers  of  the  square  root 
of  the  distance  from  the  crack  tip.  The  unknown  coefficients  in  the  equation 
were  determined  from  a  least  square  error  fit  to  the  finite  element  results. 
The  percent  error  between  the  calculated  and  the  finite  element  results  for 
compliance  was  within  ±  3%  over  the  ranges  of  0  _<  a  _<  a/W  -  0.175  and 
0.25  <_  a/W  <_  0.8.  Since  the  truncated  series  was  determined  from  curve  fits 
to  data  within  the  above  ranges,  extrapolation  outside  of  the  ranges  is  not 
recommended . 

The  second  equation  was  developed  by  adding  a  correction  term  to  an 
approximate  solution  for  the  compliance  far  from  the  crack  tip.  The  resulting 
expression  was  nonlinear  in  terms  of  the  unknown  coefficients.  An  iterative 


Crial-and-error  method  was  combined  with  a  least  square  error  method  to 
obtain  solutions  for  the  unknown  coefficients.  The  percent  error  between  the 
calculated  and  the  finite  element  results  was  within  ±  4%  over  the  ranges 
0  <  p  <  a/W  +  0.25  and  0.25  _<  a/W  <_  0.8.  Note  that  the  range  of  the  distance 
from  the  crack  tip  is  larger  than  the  range  for  the  truncated  series. 

Since  the  correction  term  was  developed  to  account  primarily  for 
compliance  in  the  vicinity  of  the  crack  tip,  the  second  equation  could  be 
applied  when  values  of  p  are  beyond  the  indicated  upper  limit,  a/W  +  0.25. 

In  an  application  where  compliance  is  measured  and  crack  length  is 
to  be  determined,  a  simple  numerical  iterative  procedure,  Newton's  method, 
can  be  applied  to  Equation  (28).  A  value  for  crack  length  can  be  readily 
determined  while  the  number  of  iterations  is  somewhat  insensitive  to  the 
initial  estimate  of  crack  length. 
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